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(57) Abstract 



A GPS reference system and corresponding method utilizing a network of reference stations and a master station (Fig. 1) Each 
reference station (102) receives GPS signals from GPS satellites (106), it observes and computes in response a pseudo-range residual for 
each GPS satellite. The master station (108) is in communication with the reference stations to receive the computed pseudo-range residuals 
and comprises a pseudo-range residual synchronizer, an ephemeris and clock correction estimator for each GPS satellite observed by the 
reference stations and a transmitter. The pseudo-range residual symAronizer computes clock differences between the reference stations so 
as to synchronize them. Each ephemeris and clock correction estimator computes an ephemeris correction and a clock correction for the 
corresponding GPS satellite. The transmitter transmits to users the computed ephemeris and clock correction. 
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WIDE AREA DIFFERENTIAL GPS REFERENCE SYSTEM AND METHOD 

This invention was supported in part by grant 93-G-004 from the Federal 
Aviation Administration. The U.S. Government may therefore have rights to 
5 this invention. 

FIELD OF INVENTION 

The present invention relates generally to differential GPS reference 
10 systems and methods. More specifically, the present invention pertains to a 
wide area differential GPS (WDGPS) system and method that employs 
common view time transfer and minimum norm solution techniques to provide 
GPS satellite ephemeris and clock corrections to GPS users. 

15 BACKGROUND OF INVENTION 

The global positioning system (GPS) is a satellite-based navigation and 
time transfer system developed by the U.S. Department of Defense. It serves 
marine, airborne and terrestrial users. 

20 The system includes GPS satellites that each transmit L1 and L2 GPS 

signals. The L1 and L2 GPS signals comprise modulated carrier signals at the 
GPS L1 and L2 frequencies (1,575 and 1,227 MHz) which are each modulated 
by a coarse-acquisition pseudo-random (C/A) code, an encrypted precision 
pseudo-random (P) code, and navigation data. The navigation data includes 

25 the transmitting GPS satellite's orbit location and clock offset from GPS time. 

Since the P code is encrypted, it is unavailable to most users. However, 
a user can determine the pseudo-range to each GPS satellite in view based on 
the time of arrival of the C/A code contained by the GPS signals transmitted 
by the GPS satellite. The user can compute a position fix based on the 

30 determined pseudo-ranges to the GPS satellites in its view and the locations 
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and clock offsets for these GPS satellites contained in the navigation data of 
the GPS signals they transmit. 

However, the U.S. Department of Defense engages in selective 
availability and deliberately misrepresents the true location and clock offsets 
of the GPS satellites in the navigation data of the GPS signals they transmit. 
Moreover, the ionosphere and the troposphere delay the arrival of these GPS 
signals at the user's location. Thus, these factors affect the accuracy and 
integrity of GPS for the user. 

In order to improve the accuracy and integrity of GPS, GPS may be 
augmented by a network of ground-based reference stations. The reference 
stations monitor the health of the GPS signals transmitted by the GPS satellites 
in their view and generate corrections to these signals which are then 
transmitted to the user. The user can accurately determine his position based 
on the GPS signals he receives from the GPS satellites and the corrections to 
these signals he receives from the reference stations. This form of GPS is 
known as differential GPS (DGPS). 

A number of systems and methods have been developed for wide area 
DGPS (WDGPS). However, all of them suffer from various problems which 
render them undesirable. 
20 One of these approaches is called position domain WDGPS and 

combines the information from multiple wide area reference stations (WRSs) 
in the position domain. Such an approach is described in "Multi-Site Real-Time 
DGPS System Using Starfix Link: Operational Results", by D. Lapucha and M. 
Huff, Proceedings of the Fifth International Technical Meeting of the Satellite 
25 Division of the Institute of Navigation, Albuquerque, Sept. 1992. 

In this approach, each WRS acts as a local area DGPS (LDGPS) 
reference station. Based on its known fixed position, it computes for each 
GPS satellite in view a corresponding pseudo-range correction to the GPS 
satellite. Each WRS then transmits a stream of LDGPS pseudo-range 
corrections to the user. The user computes pseudo-ranges to the GPS 
satellites in his view based on the GPS signals received from these GPS 
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satellites. The user then applies the pseudo-range corrections to the computed 
pseudo-ranges. If there are M WRSs, the user computes M position fixes and 
the corresponding covariance matrices. The user then uses the covariance 
matrices to form a single weighted position fix. 
5 However, this approach does not involve the estimation of any of the 

underlying states of the GPS satellites, namely GPS satellite ephemeris and 
clock errors. Nor does it involve the estimation of the clock differences 
between the different reference stations. These type of errors behave quite 
differently and could be mitigated if they were separated. Since this approach 
10 allows these errors to be combined, it unfortunately does not allow a-priori 
information about each of the underlying states to be introduced into the 
estimation process. As a result, the accuracy of each computed pseudo-range 
correction decreases rapidly with age and user to reference station separation. 
Additionally, this approach requires M independent data streams to be 
15 sent to the user if there are M different reference stations. Thus, it requires a 
rather large data bandwidth to implement. 

Another approach is called measurement domain WDGPS and combines 
the pseudo-range measurements from multiple WRSs. Such a system is 
described in "The FAA's WIB/WDGPS Testbed and Recent Test Results", by 
20 M. Lage and B. Elrod, Proceedings of the Sixth International Technical Meeting 
of the Satellite Division of the Institute of Navigation, Salt Lake City, 1993, and 
in "GPS Wide Area Augmentation System (WAAS) Testbed Results - Phase ID 
Testbed Results", by F. Haas. M. Lage, and S. Kalinowski, Proceedings of the 
Annual Meeting of the Institute of Navigation, Colorado Springs, June 1994. 
25 In this approach, each WRS computes a pseudo-range correction to 

each GPS satellite in view. These pseudo-range corrections are sent to a wide 
area master station (WMS) rather than the user. The WMS uses well known 
common view time transfer techniques to estimate the clock differences 
between the WRSs. The estimated clock differences are then removed from 
30 the pseudo-range corrections. As a result, these synchronized pseudo-range 
corrections are all referenced to a common master clock. The WMS then 
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averages for each GPS satellite the synchronized pseudo-range corrections for 
the GPS satellite. These averaged pseudo-range corrections are then 
transmitted to the user and used by the user in the same way as described 
earlier for the position domain WDGPS approach. 
5 The measurement domain WDGPS approach requires less data 

bandwidth than the position domain WDGPS approach since a single averaged 
data stream for each GPS satellite is sent to the user rather than M different 
data streams. However, in this approach as well, the pseudo-range corrections 
degrade rapidly with age and user to reference station separation since the 
10 underlying states of the GPS satellites are not estimated and therefore the user 
cannot leverage the known behavior of the underlying states to make accurate 
position fixes. This frailty is especially pronounced if the ephemeris errors 
associated with one or more of the GPS satellites is large or the ionosphere is 
introducing errors which vary spatially. 
15 Still another approach is called state space WDGPS which uses the 

information from multiple WRSs to estimate the underlying states of the GPS 
satellites, the underlying state of the ionosphere, and the underlying clock 
differences between the WRSs. Such a system is described in "Algorithms and 
Implementation of Wide Area Differential GPS", by C. Kee and B. Parkinson, 
20 Proceedings of the Fifth International Technical Meeting of the Satellite Division 
of the Institute of Navigation, Albuquerque, Sept. 1992, and also in "Wide Area 
Differential GPS (WADGPS)", by C. Kee. Stanford University, Ph.D 
Dissertation, December 1993. 

In this approach as well, each WRS computes a pseudo-range residual 
25 (i.e., correction) for each GPS satellite in view. The computed pseudo-range 
residuals are sent to a WMS which uses one large estimator to simultaneously 
compute the clock differences between the WRSs and the ephemeris and clock 
corrections for the ephemeris and clock errors of each GPS satellite in view of 
the network. 

30 The estimator in this approach must solve a system of equations which 

contains 4K + (M-1) unknowns and MK knowns, where K is the number of 
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satellites in view and M is the number of reference stations. This stems from 
the fact that each satellite has 3 unknown components of ephemeris error and 
1 unknown clock error. In addition, there are M-1 unknown offsets between the 
clocks of the WRSs clocks. At the same time, the system of equations 
5 contains MK knowns corresponding to the K pseudo-range residuals computed 
by each of the M WRSs. 

When 4K + (M-1) > MK, the system is undetermined and the WMS must 
use a minimum norm solution. When 4K + (M-1) = MK, then the system is 
exactly specified and can be readily inverted. When 4K + (M-1) < MK, then the 

10 system is over specified and the pseudo-inverse yields an estimate with 
minimum mean square error. 

Unlike the previously discussed approaches, the state space WDGPS 
approach separates the underlying states of the system, and as such, it 
enables the use of a-priori models and side information for these states. 

15 However, its complexity impedes the straight forward integration of such 
models. 

Moreover, it requires a greater data bandwidth than is necessary. This 
occurs, because the WRS network can be rather sparse and, as such, the 
system of equations used by the WMS is frequently undetermined. 

20 Specifically, rising GPS satellites that have just appeared over the horizon are 
of great concern to WDGPS users because they can greatly improve the 
geometry and accuracy of the position fix. However, the system is always 
undetermined in the case of a rising GPS satellite. Even though the minimum 
norm solution described above yields an accurate correction for a newly rising 

25 GPS satellite, the quickly varying components of the GPS satellite clock 
correction splatter into the GPS satellite ephemeris correction. Consequently, 
the ephemeris correction must be sent as frequently as the clock correction, 
thereby increasing the required data bandwidth. 

A similar approach is described in U.S. Patent No. 5,323,322 issued on 

30 June 21, 1994 to Mueller, et. al. In this approach, ephemeris and clock 
corrections for a particular GPS satellite can be provided to a user only if the 
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GPS satellite is observed by five or more WRSs. However, as indicated 
earlier, rising GPS satellites are very important to users since they provide 
added leverage for computing accurate position fixes. Thus, this approach is 
not capable of providing a user with ephemeris and clock corrections for newly 
5 rising GPS satellites. 

SUMMARY OF INVENTION 

The foregoing problems are solved by a GPS reference system and 

10 corresponding method utilizing a network of reference stations and a master 
station. Each reference station receives GPS signals from GPS satellites it 
observes and computes in response a pseudo-range residual for each GPS 
satellite it observes. The master station is in communication with the reference 
stations to receive the computed pseudo-range residuals and comprises a 

15 pseudo-range residual synchronizer, an ephemeris and clock correction 
estimator for each GPS satellite observed by the reference stations, and a 
transmitter. The pseudo-range residual synchronizer is responsive to the 
received pseudo-range residuals in order to compute clock differences between 
the reference stations and remove the clock differences from the received 

20 pseudo-range residuals so as to synchronize them. Each ephemeris and clock 
correction estimator computes an ephemeris correction and a clock correction 
for the corresponding GPS satellite in response to each synchronized pseudo- 
range residual for that GPS satellite. The transmitter transmits to users the 
computed ephemeris and clock corrections for each GPS satellite observed by 

25 the reference stations.. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a wide area differential GPS (WDGPS) system in 
30 accordance with the present invention. 
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Figure 2 provides a block diagram of a wide area reference station 
(WRS) used in the WDGPS system of Figure 1. 

Figure 3 shows the operation of the WRS of Figure 2. 

Figure 4 provides a block diagram of a wide area master station (VMS) 
5 used in the WDGPS system of Figure 1. 

Figure 5 shows the operation of the WMS of Figure 4. 

Figure 6 demonstrates the situation where a GPS satellite is observed 
by only three WRSs in the WDGPS system of Figure 1 . 

Figure 7 demonstrates the situation where a GPS satellite is observed 
1 0 by only two WRSs in the WDGPS system of Figure 1 . 

Figure 8 demonstrates the situation where a GPS satellite is observed 
by only one WRSs in the WDGPS system of Figure 1 . 

Figure 9 provides a block diagram of a user's GPS positioning unit used 
in the WDGPS system of Figure 1. 
15 Figure 10 shows the operation of the user's GPS positioning unit of 

Figure 9. 

DESCRIPTION OF SPECIFIC EMBODIMENTS 

20 Referring initially to Figure 1, there is shown a wide area differential 

GPS (WDGPS) network 100 that includes a WDGPS reference or augmentation 
system 104. The WDGPS reference system includes a network of wide area 
reference stations (WRSs) 102. Each WRS receives L1 and L2 GPS signals 
from each of the GPS satellites 106 that are in its view. In response, each 

25 WRS observes each of the GPS satellites in its view and computes a pseudo- 
range residual and an ionospheric delay estimate for each observed GPS 
satellite. The pseudo-range residuals and ionospheric delay estimates 
computed by each WRS are provided to the wide area master station (WMS) 
108 of the WDGPS reference system via a data communication link 110. The 

30 WMS in response computes ephemeris and clock corrections for each GPS 
satellite observed by the WRSs and computes a grid of ionospheric corrections. 
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The WMS then transmits the ephemeris, clock, and ionospheric corrections to 
the geostationary satellite 112. The geostationary satellite 112 re-transmits 
these corrections so that a mobile user's GPS positioning system 114 receives 
them along with the L1 GPS signals from each GPS satellite it observes. The 
5 user's GPS positioning system then computes its position based on the 
received GPS signals and the received ephemeris, clock, and ionospheric 
corrections. 

Wide Area Reference Stations (WRSs) 

10 

As shown in Figure 2, each WRS 102 includes a conventional survey 
quality dual frequency GPS receiver 116, a meteorological sensor 118, a data 
processor 120, and a data link interface 121. The GPS receiver includes a 
conventional GPS antenna 122 for receiving L1 and L2 GPS signals. The data 

15 processor includes a central processing unit (CPU) 124 that runs the pseudo- 
range smoother 126, the navigation data processor 128, the tropospheric 
correction estimator 130, the pseudo-range residual estimator 132, and the 
clock steerer 134 which are all stored in the memory 136. The memory also 
stores the precisely surveyed location 138 of the WRS. 

20 Turning to Figures 1 and 3, the GPS antenna 122 of each WRS 102 

receives L1 and L2 GPS signals from each of the GPS satellites 106 in view 
of the WRS and the GPS receiver 116 in response observes each of these 
GPS satellites. For each observed GPS satellite, the GPS receiver computes 
at each of its measurement epochs L1 and L2 pseudo-ranges p^ L1 and 

25 (from the WRS to the corresponding GPS satellite) based on the L1/CA and 
L2/CA codes received from the GPS satellite, where k„, represents the kth GPS 
satellite observed by the mth WRS. In addition, the GPS receiver makes at 
each epoch L1 and L2 carrier phase measurements 0*^ L1 and of the L1 
and L2 carrier signals transmitted by each GPS satellite it observes. In the 

30 preferred embodiment, the measurement epochs occur at one second intervals 
(i.e., 1 Hz). 

8 



WO 96/22546 



PCT/US96/00672 



The sets of L1 and L2 pseudo-ranges {pfeufc, and {p^^, and the 
sets of L1 and L2 carrier phase measurements and {0^}^, at 

each epoch are then provided to the pseudo-range smoother 122 of the WRS, 
where K„ represents the total number of GPS satellites observed by the mth 
5 WRS. For each GPS satellite observed by the WRS, the pseudo-range 
smoother uses the corresponding L1 and L2 earner phase measurements and 
the corresponding L1 and L2 pseudo-ranges to produce at each epoch a 
corresponding ionospheric delay estimate ft and a corresponding smoothed 
pseudo-range p*™ from which the ionospheric delay estimate has been 

10 removed. This is done using the method described in "GPS Survey", by C. 
Goad, a chapter in The Global Positioning System: Theory and Applications", 
edited by B. Parkinson, J Spilker, P. Axelrod, and P. Enge, American Institute 
of Aeronautics and Astronautics, 1995, which is hereby explicitly incorporated 
by reference. The set of smoothed pseudo-ranges {p^^i at each epoch are 

15 sent to the pseudo-range residual estimator 132 and the set of ionospheric 
delay estimates {O 1 ^,, at each epoch are sent to the WMS 108. 

The GPS receiver 116 of each WRS 102 also extracts the navigation 
data contained by the L1 and L2 GPS signals it receives from each GPS 
satellite it observes. The navigation data processor 128 is provided with the 

20 navigation data for each observed GPS satellite and unpacks from it the GPS 
satellite's location vector r** B and the difference (or offset) B"™ 6 from GPS time 
of tne GPS satellite's clock. Then, for each GPS satellite observed by the 
WRS, the navigation data processor computes at each epoch a line of sight 
vector 1^and a nominal range r*™- 8 • l^from the WRS to the GPS satellite 

25 based on the unpacked GPS satellite's location and the known reference 
station location 138 stored in the memory 130. The set of line of sights 
{I'Jlfci are provided to the WMS 108 while the set of computed nominal 
ranges {r*- B • l 1 *} 1 ^, and the set of unpacked clock differences {B fcw,B }'!* Bl are 
provided to the pseudo-range residual estimator 132. 

30 The conventional meteorological sensor 118 of each WRS 102 senses 

surface pressure, humidity, and/or temperature and provides this sensed 

9 
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meteorological data to the tropospheric correction estimator 130. The 
tropospheric delay estimator 130 periodically computes at each epoch a 
tropospheric delay estimate T^for each of the GPS satellites 106 observed by 
the WRS in response to the sensed meteorological data and the line of sight 
5 to each GPS satellite. This is done using a modified Hopfield method as 
described in "Analysis of Tropospheric Delay Prediction Models: Comparisons 
with Ray-Tracing and Implications for GPS Relative Positioning", by H. Janes, 
R Langley, and S. Newby, Bulletin Geodesique, Vol 65, pp 151-161, 1991, 
which is hereby incorporated by reference. The set of tropospheric delay 
10 estimates {T^ l ^ 1 are then provided to the pseudo-range residual estimator 
132. 

Each smoothed pseudo-range received at each epoch by the 
pseudo-range residual estimator 132 may be modeled as: 

Pm=r k -1^ + b m -B^ + T m km + vJ: Eq.(1) 

where r*™ is a vector that represents the actual location of the kth GPS satellite 
15 106 observed by the WRS, b m represents the actual difference from GPS time 
of the clock of the GPS receiver 116 of the mth WRS 102, represents the 
actual difference from GPS time of the clock of the kth GPS satellite observed 
by the mth WRS, and v^ represents the noise associated with the kth GPS 
satellite and the mth WRS that observes it. For each smoothed pseudo-range 
20 p^, the pseudo-range residual estimator 132 subtracts the corresponding 
nominal range r*™* • 1^, unpacked clock offset B 1 ^ 8 , and tropospheric 
correction T^to produce a pseudo-range residual (i.e., correction) Ap^at each 
epoch for the error in the smoothed pseudo-range. This pseudo-range residual 
may be represented as: 

Ap^T = Ar k -lSb m - AB^ + vS: *M2> 

25 where Ar*™ = r 1 ^ - i* 1 * 3 and is a vector representing the correction for the 
ephemeris error of the kth GPS satellite observed by the mth WRS, and AB^ 

10 
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■ - B lwB and represents the correction for the clock error of the kth GPS 
satellite observed by the mth WRS. 

Since the clock difference b m dominates the expression in Equation (2), 
the reference station clock steerer 134 averages at each epoch all of the 
5 pseudo-range residuals {Ap^^., for the GPS satellites observed by the WRS 
to form a rough estimate b m of the actual clock difference b m , as is well known 
by those skilled in the art. The clock difference estimate b m j S subtracted by 
the pseudo-range residual estimator 132 from each pseudo-range residual A ft *" 
to produce an adjusted pseudo-range residual Ap^ given by: 

Ap^Ar'-l* +Ab m -AB k + v ;, Eq.(3) 

10 where Ab m ■ b m - b m . The set of pseudo-range residuals {Ap^ 1 ^ are then 
transmitted to the WMS 108 over one of the data communication links 1 1 0 by 
the data link interface 121. Each data communication link 110 may comprise 
a conventional modem/telephone line system, an RF or other frequency 
transmitter/receiver system, a fiber optic cable network, or other suitable 

15 transmission system. 

Wide Area Master Station (WMS) 

Referring to Figure 4, the WMS 108 includes data link interfaces 140, a 
20 data processor 142, and a transmitter 143. The data processor includes a 
pseudo-range residual synchronizer 144, ephemeris and clock correction 
estimators 146, and an ionospheric grid generator 148 which are run by the 
CPU 147 and stored in the memory 150. 

Turning to Figures 1 and 5, the WMS 108 includes the same number of 
25 data link interfaces 140 as there are WRSs 102 in the WDGPS reference 
system 104. As a result, the interfaces 140 provide the pseudo-range residual 
synchronizer 144 at each epoch with the pseudo-range residuals {{Ap^^JJJ,,, 
provide the ionospheric grid generator 148 with the ionospheric delay estimates 
{{I'S'&.iJm.i. and provide the ephemeris and clock correction estimators 146 
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with the lines of sight {{1^)1, (re-phrased as {{iSjE-ift-i. as wi » be 
discussed later), where M represents the total number of WRSs. 

As those skilled in the art will appreciate, the lines of sight {{1^^™.,}"., 
could be computed by the WMS 108 rather than the WRSs 102. In this case, 
5 the WMS may include a GPS receiver that receives the GPS signals 
transmitted by the GPS satellites 104 and then removes the navigation data 
contained in the GPS signals. Or, the navigation data in the GPS signals 
received by the WRSs could be simply provided to the WMS via the 
communication links 110. In either case, based on the GPS satellite locations 

10 contained in the navigation data of the GPS signals and the known locations 
of the WRSs which would be stored in the memory 150 of the WMS, the WMS 
would then compute the lines of sight. 

In response to the ionospheric delay estimates {{I^X,}"*,, the 
ionospheric grid generator 148 produces a grid of ionospheric corrections 

15 {IJ^T over the area covered by the GPS satellites 106. This is done in the 
manner described in "Flight Trials of the Wide-Area Augmentation System 
(WAAS)", by Walter et. al., Proceedings of the International Technical Meeting 
of the Satellite Division of the Institute of Navigation, Salt Lake City, September 
1994, which is hereby explicitly incorporated by reference. The grid of 

20 ionospheric corrections is provided to the user's positioning system 1 14 in a 
correction message transmitted by the transmitter 143. 

In response to the pseudo-range residuals {{Ap^,}".,, the pseudo- 
range synchronizer 144 first performs at each epoch a common view time 
transfer to estimate the differences between the clocks of the WRSs 102 with 

25 respect to a master clock. The clock differences are then removed from the 
pseudo-range residuals so as to synchronize them. 

In order to do so, one of the WRSs is designated as a primary WRS and 
the remaining WRSs are designated as secondary WRSs. For each GPS 
satellite 106 in common view of the primary WRS and at least one of the 

30 secondary WRSs, the pseudo-range residual synchronizer 144 first subtracts 
at each epoch the pseudo-range residual Apfr M computed for it by the primary 
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WRS from each pseudo-range residual A~w 

PnT u computed for it by those of the 

secondary WRSs that observe It, where M represents the primary WRS (in 

addition to representing the number of WRSs) and. for m * M, represents 

the kth of the GPS satellites that is observed by the mth secondary WRS and 

5 by the primary WRS. This results in the difference A^"» at each epoch which 

is given as follows: 

Eq. (4) 

aSt = Ajj- - ApJr - Ar*~ . c£- - ijn ♦ Ab m - Ab M ♦ -yssi 

where vjj-jj = - ift* Assuming that the difference in line of sights 
OSF* - 1m M ) «s small and/or the e'phemeris correction Ar**" is small, then, an 
estimate of the WRS clock error can be made by averaging the difference A^ M 
10 as follows: 

Abm " AbM 85 k^; 2 ^ ^ Eq - (5) 

where Ab m - Ab M represents the clock difference at each epoch between the 
primary WRS and the mth secondary WRS and !<„, „ represents the total 
number of GPS satellites in common view of the primary WRS and the mth 
WRS. Moreover, if the clocks of the GPS receivers 1 1 6 of the WRSs are high 
15 quality clocks, such as rubidium clocks, then the difference Ab m - Ab M can be 
averaged over time to take advantage of the well behaved drift of the clocks. 
In this case, the expression in Equation (5) can be re-written as: 

Ab m - Ab M . JL h ty , A ^ M (t|) Eq (6) 

where h(tj) represents the well known impulse response for a low pass filter. 
From the foregoing, it is evident that in order to compute a clock 
20 difference Ab m - ADm at each epoch for each WRS 102, the secondary WRSs 
must be positioned with respect to the primary WRS so that each secondary 
WRS will always have at least one GPS satellite in common view with the 
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primary WRS. Preferably, the WRS that is closest to the WMS and most 
centrally located with respect to the more peripheral WRSs is the primary WRS. 

The pseudo-range synchronizer 144 then at each epoch subtracts from 
the pseudo-range residuals {Ap*^, computed by each of the secondary 
5 WRSs the corresponding clock difference Ab m - Ab M to produce the pseudo- 
range residuals {Apftfa, at each epoch. Thus, the pseudo-range residuals 
computed by the secondary WRSs {{Apftfc-id are synchronized with the 
pseudo-range residuals {Apfc} 1 ^ computed by the primary WRS and 
referenced to the same master clock (i.e., the clock of the primary WRS). As 
10 a result, each synchronized pseudo-range residual can now be represented as 
Ap^and given by: 

A? m m = Ar"" • ft ♦ AB^ + Eq.(7) 

for m = 1 to M, and where Apfr = Apft and AB*™ = Ab M - AB^ and represents 
the clock correction for the kth GPS satellite observed by the mth WRS 
synchronized to the clock of the primary WRS. 
15 Furthermore, Equation (7) can be re-written as follows: 

A? mh = Ar k -1^ + AB k + v^ Eq.(8) 

where k represents the kth of all of the GPS satellites 106 observed by the 
entire network of WRSs 102 and m„ represents the mth of those of the WRSs 
that observe the kth GPS satellite. The set of synchronized pseudo-range 
residuals {{Ap^,}^,, and the corresponding set of lines of sight {{1JXi)£.i 
20 (re-phrased to correspond with the synchronized pseudo-range residuals) are 
then provided to the ephemeris and clock correction estimators 146 of the 
WMS 108, where M k represents the total number of WRS that observe the kth 
GPS satellite. 

Since the clock differences of the WRSs 102 have been removed from 
25 the set of synchronized pseudo-range residuals {{ApHjE.,}^,, the WMS 108 
can employ an ephemeris and clock correction estimator 146 for each of the 
GPS satellites 106 observed by the WRSs. Each ephemeris and clock 
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correction estimator then computes at each epoch the ephemeris correction 
vector Ar* and the clock correction AB k for the corresponding GPS satellite in 
response to the lines of sight {1jy^, to the corresponding GPS satellite and 
the synchronized pseudo-range residuals {Ap^J^ 1 for the corresponding GPS 
5 satellite, independent of the synchronized pseudo-range residuals for the other 
GPS satellites. As will be evident from the following discussion, by employing 
separate ephemeris and clock correction estimators 146, the computational 
complexity of a single large estimator is avoided. 

10 Four Dimensional Estimators 

In one embodiment, each ephemeris and clock correction estimator 146 
solves for the four dimensional state vector [Ar*, AB k ] T based on the following 
matrix system: 



r . -k i 
Ap lk 










. -y 




(0 T -1 


Ar k 
AB k 


Eq. (9) 













1 5 The system of Equation (9) can be more compactly re-written as follows: 

Ap k =G w x k Eq.(10) 

where Ap k is the vector containing each synchronized pseudo-range residual 
Ap* computed by the WRSs 102 for the kth GPS satellite 106, x k is the four 
dimensional state vector containing the ephemeris and clock corrections Ar* 
and AB k for the kth GPS satellite, and Gw is the design matrix that contains, for 
20 each WRS that observes the kth GPS satellite, a line of sight from it to the kth 
GPS satellite and a matrix element (i.e., -1) that indicates that the 
corresponding synchronized pseudo-range residual in the vector Ap k should be 
adjusted by the clock correction AB k for the kth GPS satellite. 
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When the WDGPS reference system 104 includes at least four WRSs 
102 (i.e., the primary WRS and at least three secondary WRSs) and the kth 
GPS satellite 106 is observed by exactly four of these WRSs t then the system 
of Equation (10) is exactly specified. Furthermore, when the WDGPS reference 
5 system 104 includes more than four WRSs 102 and the kth GPS satellite 106 
is observed by more than four of these WRSs, then the system of Equation 
(10) is over specified. In either case, the corresponding estimator 146 first 
computes a pseudo-inverse {G^G^^G^ of the design matrix In the case 
where the corresponding GPS satellite is observed by exactly four of the 
1 0 WRSs, the pseudo-inverse reduces to the inverse G^ of the design matrix G^ 
The estimator then computes a minimum mean square error solution for the 
state vector x k as follows: 

x k = (Gw G w r 1 Gyl Ap k Eq.(11) 

Moreover, the system of Equation (10) is under specified when the kth 
GPS satellite 106 is observed by less than four WRSs 102, regardless of 
1 5 whether the WDGPS reference system 1 04 includes less than four WRSs (i.e., 
the primary WRS and no more than two secondary WRSs) or at least four 
WRSs. In this situation, the corresponding estimator 146 computes a pseudo- 
inverse (GvjGwV'Gj of the design matrix G w and then computes a minimum 
norm solution for the state vector x k as follows: 

x k = Gw (G w Gw)' 1 Ap k Eq. (12) 

20 The minimum norm solution has the minimum length or energy among all of 
the satisfactory solutions. 

The set of computed ephemeris and clock corrections {Ai*; AB 1 *}^ for 
the GPS satellites observed by the WRSs 102 are then provided to the 
transmitter 144. The transmitter transmits them as corrections to the user's 

25 GPS positioning system 1 14 via the geostationary satellite 1 12 or other suitable 
broadcast facility so that the user's positioning system can compute an 
accurate position fix. Alternatively, as those skilled in the art will appreciate, 
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the computed ephemeris and clock corrections may be transmitted directly to 
the user's positioning system by the transmitter 143. 

Figures 6 illustrates the situations when the kth GPS satellite 106 is 
observed by only three WRSs 102 and the system of Equation (10) is under 
5 specified. Since the kth GPS satellite 106 is observed by only three WRSs, 
this geometry does not provide leverage for precisely observing each of the 
radial (i.e., orbit radius), along track (along satellite velocity), and cross track 
(i.e., orthogonal to radial and along track) components of the ephemeris error 
of the GPS satellite and the clock error of the GPS satellite. However, when 

10 the user's positioning system 1 14 is between the three WRSs or adjacent to 
one of the three WRSs, the minimum norm solution of Equation (12) does 
provide accurate ephemeris and clock corrections Ar* and AB* for the 
ephemeris and clock errors since the components of these errors that effect the 
user's positioning system effect the three WRSs substantially in the same 

15 manner. 

Figure 7 demonstrates the case when the kth GPS satellite 106 is 
observed by only two WRSs 102 and the system of Equation (10) is under 
specified. Similar to the situation described for Figure 6, although this 
geometry does not provide leverage for precisely observing each of the 

20 components of the ephemeris and clock errors, the minimum norm solution of 
Equation (12) does provide accurate ephemeris and clock corrections Ar* and 
AB k for when the user's positioning system 114 is between the two WRSs or 
adjacent to one of the WRSs. As in the situation of Figure 6, this is due to the 
fact that the components of the ephemeris and dock errors that effect the 

25 user's positioning system effect the two WRSs substantially in the same 
manner. 

Similarly, Figure 8 shows the situation when the kth GPS satellite 106 
is observed by only one WRS 102 and the system of Equation (10) is under 
specified. This geometry also does not provide the leverage for accurately 
30 observing the components of the ephemeris and clock errors. But, when the 
user's positioning system 114 is adjacent to the WRS, the minimum norm 
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solution of Equation (12) provides accurate ephemeris and clock corrections 
Ar* and AB k because the components of the ephemeris and clock errors that 
effect the user's positioning system effect the WRS in substantially the same 
manner. 

5 From the foregoing, it is evident that the WDGPS reference system 104 

is capable of estimating applicable ephemeris and clock errors for any GPS 
satellite in view of at least one of the WRSs 102, provided that the user's 
positioning system 114 is between WRSs or adjacent to a WRSs. As result, 
the WDGPS reference system can provide to users accurate and valid 
10 ephemeris and clock corrections for each GPS satellite from the time that it is 
first observed by one of the WRSs (i.e., appears on the horizon) until the time 
that it is no longer observed by any of the WRSs (i.e., disappears below the 
horizon). 

However, in this embodiment, even though the four dimensional 
15 minimum norm solution for the state vector x k always generates corrections 
which are valid for users in between or adjacent the WRSs 102, the high 
frequency energy in the clock correction AB k (resulting from the high frequency 
energy in the clock error due to selective availability) splatters into the 
ephemeris corrections Ar\ As a consequence, the ephemeris corrections 
20 {Ar*}^, must be transmitted by the WMS 108 in the correction message to 
users as frequently (i.e., at each epoch) as the clock corrections {AB k }" Bl . This 
may be undesirable in an operational WDGPS reference system where the 
broadcast bandwidth is severely limited and only allows the clock corrections 
to be sent frequently. In such a case, three dimensional ephemeris and clock 
25 correction estimators may be used instead. 
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Three Dimensional Estimators 

In this embodiment, each ephemeris and clock correction estimator 146 
first subtracts each synchronized pseudo-range residual ApJ^ for the 
5 corresponding GPS satellite computed by the secondary WRSs from the 
synchronized pseudo-range residual Ap^ for the corresponding GPS satellite 
computed by the primary WRS, where M k also represents the primary WRS (in 
addition to representing the total number of WRSs that observe the kth GPS 
satellite). Each difference Aj^ can be expressed as: 

Eq. (13) 

AlU = a p5» " *A = Ar^ • (1^ - ijy ♦ 

10 where = vj^ - 1^. This is done in order to separate the computation of 
the clock correction AB k from the computation of the ephemeris correction Ar* 
so that the high frequency energy of the clock correction does not splatter into 
the ephemeris correction. 

Each estimator then solves the three dimensional ephemeris correction 

15 vector Ar" based on the following matrix system: 



r A k 




d!. 


- it/ 


A* 

m n» M n 








A k 




0t,-i 


- V 



Similar to the system of Equation (9), the system of Equation (14) can be more 
compactly re-written as: 

A k = G w Ar k Eq.(15) 

where A k is a vector containing each computed difference Aj^ and G w is the 
matrix that contains, for each secondary WRS that observes the kth GPS 
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satellite, the difference between the line of sight 1^ from the corresponding 
secondary reference station to the kth GPS satellite and the line of sight 1^ 
from the primary reference station to the kth GPS satellite. 

As with the system of Equation (10), the system of Equation (15) is 
5 exactly specified when the WDGPS reference system 104 includes at least four 
WRSs 102 (i.e., the primary WRS and at least three secondary WRSs) and the 
kth GPS satellite 106 is observed by exactly four of these WRSs. And, the 
system of Equation (15) is over specified when the WDGPS reference system 
104 includes more than four WRSs 102 and when the kth GPS satellite 106 is 
10 observed by more than four of these WRSs. The corresponding estimator 146 
in this case computes a three dimensional minimum mean square error solution 
for the state vector Ar* as follows: 

Ar k = (GwGw)" 1 Gw A k Eq.(16) 

In addition, the system of Equation (15) is under specified if the kth GPS 
satellite 106 is observed by less than four WRSs 102 in either situation where 
15 the WDGPS reference system 104 includes less than four WRS (i.e., the 
primary WRS and no more than two secondary WRSs) or at least four WRSs. 
In this case, the corresponding estimator 146 computes a three dimensional 
minimum norm solution for the state vector Ar" as follows: 

Ar k = Gyl (G w GyJ) 1 A k Eq-(17) 

For each of the situations described above, after computing the 
20 ephemeris correction Ar*, each estimator 146 then computes the clock 
correction AB k according to: 

The set of ephemeris and clock corrections {Ar*; AB k }" Bl computed by 
the estimators 146 are then provided to the transmitter 144 for transmission as 
ephemeris and clock corrections. 
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As was alluded to earlier, the three dimensional estimators 146 prevent 
the high frequency energy in the clock correction AB k from splattering into the 
computed ephemeris corrections Ar\ As a result, the set of computed 
ephemeris corrections {Ar*}^ need not be transmitted as often as the set of 
5 computed clock corrections 

However, since newly rising GPS satellites are important to GPS users, 
a disadvantage to this three dimensional approach is that corrections for newly 
rising GPS satellites cannot be computed until the GPS satellite is viewed by 
the primary WRS and at least one secondary WRS. For this reason, the 
10 estimators 146 may utilize the three dimensional approach for each GPS 
satellite that is observed by the primary WRS and at least one secondary WRS 
and utilize the four dimensional approach for each GPS satellite that is 
observed by at least one WRS but is not observed by the primary WRS and at 
least one secondary WRS. 

15 

User's GPS Positioning Unit 

Referring to Figure 1, since the WDGPS reference system 104 is 
capable of estimating the ephemeris and clock corrections for any GPS satellite 

20 in view of at least one of the WRSs 102. the user's GPS positioning unit 114 
can make accurate position fixes using the ephemeris and clock corrections 
transmitted by the WMS 108. 

As shown in Figure 9, the user's positioning unit 114 includes a 
correction message receiver 152, a conventional L1 frequency GPS receiver 

25 154, a data processor 158, and a display 160. The data processor includes 
a CPU 166 that runs the correction message unpacker 168. the pseudo-range 
correction estimator 170, the pseudo-range estimator 172, the ionospheric 
correction estimator 174, the tropospheric correction estimator 176, and the 
position estimator 178 which are all stored in the memory 180. The receiver 

30 1 52 includes an antenna 162 for receiving the correction message transmitted 
by the WMS 108 and re-transmitted by the geostationary satellite 112. The 
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GPS receiver 154 includes a conventional GPS antenna 164 for receiving L1 
GPS signals from GPS satellites in its view. Alternatively, the correction 
message receiver 152 and the GPS receiver 154 may be integrated as a single 
unit for receiving both the correction message and the GPS signals. 
5 Turning to Figures 1 and 10, the GPS antenna 164 of the GPS receiver 

154 receives an L1 GPS signal from each of the GPS satellites 106 that is in 
view and the GPS receiver 154 in response observes each of these GPS 
satellites. For each observed GPS satellite, the GPS receiver 154 computes 
at each epoch a pseudo-range from the user's positioning unit 1 14 to the 

1 0 corresponding GPS satellite based on the L1/CA code received from the GPS 
satellite, where k„ represents the kth GPS satellite observed by the user's 
positioning unit. At each epoch, the GPS receiver 154 also computes a line of 
sight 1J^ from the user's positioning unit to each observed GPS satellite. The 
set of pseudo-ranges {p^, at each epoch are then provided to the pseudo- 

15 range estimator 172 of the user's positioning unit and the set of lines of sight 
{ 1 2"&i at each epoch are provided to the tropospheric correction estimator 
176. where K„ represents the total number of GPS satellites observed by the 
user's positioning unit. 

The tropospheric correction estimator 176 computes a tropospheric 

20 correction T^-for each of the GPS satellites 106 observed by the user's 
positioning unit in response to the line of sight 1j» to each GPS satellite. This 
is done according to the modified Hopfield method discussed earlier. The set 
of tropospheric corrections {T^}^ are then provided to the pseudo-range 
estimator 172. 

25 Alternatively, the user's positioning system could include a 

meteorological sensor to sense surface pressure, humidity, and/or temperature. 

In response to this meteorological data, the tropospheric correction estimator 

176 would compute the tropospheric corrections {Tj| u }|^, 1 . 

The antenna 162 of the receiver 152 receives the correction message 
30 transmitted by the WMS 108 and the receiver in response provides the 

correction message to the correction message unpacker 168. The correction 
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message unpacker unpacks the sets of ephemeris and clock corrections 
{Ar"; AB*}^ and sends them to the pseudo-range estimator 172. In addition, 
it unpacks the grid of ionospheric corrections {l$*£ n and sends them to the 
ionospheric delay estimator 174. 
5 In response to the grid of ionospheric corrections {IJ'tT. the ionospheric 

delay estimator computes an ionospheric correction for each GPS satellite 
observed by the user's positioning unit. The set of ionospheric corrections 
ar© provided to the pseudo-range estimator 172. 
Alternatively, the user's positioning system could use a dual frequency 
10 GPS receiver to receive L1 and 12 GPS signals. In this case, the dual 
frequency GPS receiver could compute the ionospheric corrections directly 
without requiring the grid of ionospheric corrections transmitted by the WMS 
108 to compute these corrections. 

In response to the unpacked ephemeris and clock corrections Ar*" 
1 5 and AB"" for each GPS satellite observed by the user's positioning system, the 
pseudo-range correction estimator 172 computes at each epoch a pseudo- 
range correction Ap^ for the corresponding GPS satellite according to: 

ApS" = a£ • £ - AB"- Eq.(19) 

and computes at each epoch the term (t - 1,) Ap£\ where A/£ is the pseudo- 
20 range correction velocity for the corresponding GPS satellite computed by 
differencing the computed pseudo-range corrections from the last two epochs 
as follows: 

A .ju = Apfrg - Apfo.,) Eq> (20) 

Each computed pseudo-range correction Apj^ , term (t - Ap k , 
tropospheric delay estimate TjJ", and ionospheric correction is subtracted 
25 from the corresponding pseudo-range p^by the pseudo-range estimator 1 72 to 
produce the corrected pseudo-range /ftgiven by: 
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Pu - Pu ' ApJ* - (t - t,) Ap^ - 1^ - T u *" Eq- (2D 

The set of corrected pseudo-ranges {/ft}!^ are then provided to the 
position estimator 178. The position estimator computes in response a position 
fix for the user's positioning unit 1 14 according to well known techniques. The 
computed position fix is provided to the display 160 which then displays the 
5 position fix for the user. 

While the present invention has been described with reference to a few 
specific embodiments, the description is illustrative of the invention and is not 
to be construed as limiting the invention. Various modifications may occur to 
those skilled in the art without departing from the true spirit and scope of the 
10 invention as defined by the appended claims. 
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What is claimed is: 

1 . A GPS reference system for providing GPS satellite ephemeris and clock 
corrections to users, the GPS reference system comprising: 

5 a network of reference stations, each reference station receiving GPS 

signals from GPS satellites it observes and computing in response a pseudo- 
range residual for each GPS satellite rt observes; and 

a master station in communication with the reference stations to receive 
the computed pseudo-range residuals, the master station comprising: 
10 a pseudo-range residual synchronizer, responsive to the received 

pseudo-range residuals, to compute clock differences between the reference 
stations and remove the clock differences from the received pseudo-range 
residuals so as to synchronize them in time; 

for each GPS satellite observed by the reference stations, a 
15 corresponding ephemeris and clock correction estimator, responsive to each 
synchronized pseudo-range residual for the corresponding GPS satellite, to 
compute an ephemeris correction and a clock correction for the corresponding 
GPS satellite; and 

a transmitter to transmit to users the computed ephemeris and 
20 clock corrections for each GPS satellite observed by the reference stations. 

2. A GPS reference system as recited in claim 1 wherein each ephemeris 
and clock correction estimator is configured to compute the ephemeris and 
clock corrections for the corresponding GPS satellite when the corresponding 

25 GPS satellite is observed by at least one of the reference stations. 

3. A GPS reference system as recited in claim 2 wherein the network of 
reference stations comprises no more than three reference stations. 

30 4. A GPS reference system as recited in claim 2 wherein the network of 
reference stations comprises at least four reference stations. 
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5. A GPS reference system as recited in claim 1 wherein: 

the network of reference stations comprises no more than three 
reference stations; 

each ephemeris and clock correction estimator is configured to compute 
5 the ephemeris and clock corrections for the corresponding GPS satellite in 
accordance with: 

x* = Gj(G w G v Ir 1 Ap k 

where: 

k represents the corresponding GPS satellite; 

x* is a state vector containing the ephemeris and clock corrections for 

10 the corresponding GPS satellite; 

Ap k is a vector containing each synchronized pseudo-range residual for 
the corresponding GPS satellite, each synchronized pseudo-range residual in 
the vector Ap k corresponding to one of the reference stations that observes the 
corresponding GPS satellite; and 

15 Gw is a matrix containing, for each reference station that observes the 

corresponding GPS satellite, a line of sight from it to the corresponding GPS 
satellite and an element that indicates that the corresponding synchronized 
pseudo-range residual in the vector Ap k should be adjusted by the clock 
correction for the corresponding GPS satellite. 

20 

6. A GPS reference system as recited in claim 1 wherein: 

the network of reference stations comprises no more than three 
reference stations and includes a primary reference station and no more than 
two secondary reference stations; 
25 each ephemeris and clock correction estimator is configured to compute 

the ephemeris correction for the corresponding GPS satellite when it is 
observed by the primary reference station and at least one of the secondary 
reference stations in accordance with: 
where: 
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Ar' = Gvf (G w GvJr 1 A* 

k represents the corresponding GPS satellite; 
Ar* is a vector containing the ephemeris correction for the corresponding 
GPS satellite; 

A k is a vector containing, for each secondary reference station that 
5 observes the corresponding GPS satellite, the difference between the 
synchronized pseudo-range residual for the corresponding GPS satellite 
computed by it and the synchronized pseudo-range residual for the 
corresponding GPS satellite computed by the primary reference station; 

Gw is a matrix containing, for each secondary reference station that 
10 observes the corresponding GPS, the difference between a line of sight from 
it to the corresponding GPS satellite and a line sight from the primary reference 
station to the corresponding GPS satellite; and 

each ephemeris and clock correction estimator is further configured to 
compute the clock correction for the corresponding GPS satellite in response 
15 to the computed ephemeris correction for the corresponding GPS satellite. 

7. A GPS reference system as recited in claim 1 wherein: 

the network of reference stations comprises at least four reference 
stations; 

20 each ephemeris and clock correction estimator is configured to compute 

the ephemeris and clock corrections for the corresponding GPS satellite when 
it is observed by no more than three of the reference stations In accordance 
with: 

and compute the ephemeris and clock corrections for the corresponding GPS 
25 satellite when it is observed by at least four of the reference stations in 
accordance with: 
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where: 

k represents the corresponding GPS satellite; 

x* is a state vector containing the ephemeris and clock corrections for 
the corresponding GPS satellite; 
5 Ap k is a vector containing each synchronized pseudo-range residual for 

the corresponding GPS satellite, each synchronized pseudo-range residual in 
the vector Ap k corresponding to one of the reference stations that observes the 
corresponding GPS satellite; 

Gyy is a matrix containing, for each reference station that observes the 
10 corresponding GPS satellite, a line of sight from it to the corresponding GPS 
satellite and an element that indicates that the corresponding synchronized 
pseudo-range residual in the vector Ap k should be adjusted by the clock 
correction for the corresponding GPS satellite. 

15 8. A GPS reference system as recited in claim 1 wherein: 

the network of reference stations comprises at least four reference 
stations including a primary reference station and at least three secondary 
reference stations; 

each ephemeris and clock correction estimator is configured to compute 
20 the ephemeris correction for the corresponding GPS satellite when it is 
observed by no more than three of the reference stations in accordance with: 

Ar k =G v I(G w G v Ir 1 A* 

and compute the ephemeris correction for the corresponding GPS satellite 
when it is observed by at least four of the reference stations in accordance 
with: 

Ar k = (G^G*)' 1 Gj A* 

25 where: 
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k represents the corresponding GPS satellite; 
Ar* is a vector containing the ephemeris correction for the corresponding 
GPS satellite; 

A* is a vector containing, for each secondary reference station that 
5 observes the corresponding GPS satellite, the difference between the 
synchronized pseudo-range residual for the corresponding GPS satellite 
computed by it and the synchronized pseudo-range residual for the 
corresponding GPS satellite computed by the primary reference station; and 
G w is a matrix containing, for each secondary reference station that 
10 observes the corresponding GPS, the difference between a line of sight from 
it to the corresponding GPS satellite and a line sight from the primary reference 
station to the corresponding GPS satellite; and 

each ephemeris and clock correction estimator is further configured to 
compute the clock correction for the corresponding GPS satellite in response 
15 to the computed ephemeris correction for the corresponding GPS satellite. 

9. A GPS reference system for providing GPS satellite ephemeris and clock 
corrections to users, the GPS reference system comprising: 

a network of reference stations including a primary reference station and 

20 at least one secondary reference station, each reference station receiving GPS 
signals from GPS satellites it observes and computing in response a pseudo- 
range residual for each GPS satellite it observes, each secondary reference 
station positioned with respect to the primary reference station such that at 
least one of the GPS satellites it observes is simultaneously observed by the 

25 primary reference station; and 

a master station in communication with the reference stations to receive 
the computed pseudo-range residuals, the master station comprising: 

a pseudo-range residual synchronizer to compute clock 
differences between the reference stations in response to the received 

30 pseudo-range residuals for each GPS satellite that is simultaneously observed 
by the primary reference station and at least one of the secondary reference 
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stations and to remove the clock differences from the received pseudo-range 
residuals so as to synchronize them; 

for each GPS satellite observed by the reference stations, a 
corresponding ephemeris and clock correction estimator, responsive to each 
5 synchronized pseudo-range residual for the corresponding GPS satellite, to 
compute an ephemeris correction and a clock correction for the corresponding 
GPS satellite; and 

a transmitter to transmit to users the computed ephemeris and 
clock corrections for each GPS satellite observed by the reference stations. 

10 

10. A GPS reference system as recited in claim 9 wherein each ephemeris 
and clock correction estimator is configured to compute the ephemeris and 
clock corrections for the corresponding GPS satellite when the corresponding 
GPS satellite is observed by at least one of the reference stations. 

15 

11. A GPS reference system as recited in claim 10 wherein the network of 
reference stations comprises no more than three reference stations. 

12. A GPS reference system as recited in claim 10 wherein the network of 
20 reference stations comprises at least four reference stations. 

13. A GPS reference system as recited in claim 9 wherein: 

the network of reference stations comprises no more than three 
reference stations; 

25 each ephemeris and clock correction estimator is configured to compute 

the ephemeris and clock corrections for the corresponding GPS satellite in 
accordance with: 
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where: 

k represents the corresponding GPS satellite; 

x* is a state vector containing the ephemeris and clock corrections for 
the corresponding GPS satellite; 
5 Ap k is a vector containing each synchronized pseudo-range residual for 

the corresponding GPS satellite, each synchronized pseudo-range residual in 
the vector Ap k corresponding to one of the reference stations that observes the 
corresponding GPS satellite; and 

Gw is a matrix containing, for each reference station that observes the 
10 corresponding GPS satellite, a line of sight from it to the corresponding GPS 
satellite and an element that indicates that the corresponding synchronized 
pseudo-range residual in the vector Ap K should be adjusted by the clock 
correction for the corresponding GPS satellite. 

15 14. A GPS reference system as recited in claim 9 wherein: 

the network of reference stations comprises no more than three 
reference stations; 

each reference station is further configured to compute a line of sight to 
each of the GPS satellites it observes; 
20 the master station is also in communication with the reference stations 

to receive the computed lines of sights; 

each ephemeris and clock correction estimator is configured to compute 
the ephemeris correction for the corresponding GPS satellite when it is 
observed by the primary reference station and at least one of the secondary 
25 reference stations in accordance with: 

Ar k = Gj (GwGJ)- 1 A* 

where: 

k represents the corresponding GPS satellite; 
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Ar* is a vector containing the ephemeris correction for the corresponding 
GPS satellite; 

A k is a vector containing, for each secondary reference station that 
observes the corresponding GPS satellite, the difference between the 
5 synchronized pseudo-range residual for the corresponding GPS satellite 
computed by it and the synchronized pseudo-range residual for the 
corresponding GPS satellite computed by the primary reference station; 

Gw is a matrix containing, for each secondary reference station that 
observes the corresponding GPS, the difference between a line of sight from 
10 it to the corresponding GPS satellite and a line sight from the primary reference 
station to the corresponding GPS satellite; and 

each ephemeris and clock correction estimator is further configured to 
compute the clock correction for the corresponding GPS satellite in response 
to the computed ephemeris correction for the corresponding GPS satellite. 

15 

15. A GPS reference system as recited in claim 9 wherein: 

the network of reference stations comprises at least four reference 
stations; 

each ephemeris and clock correction estimator is configured to compute 
20 the ephemeris and clock corrections for the corresponding GPS satellite when 
it is observed by no more than three of the reference stations in accordance 
with: 

x k = GvJ (G w GvJ)- 1 Ap k 

and compute the ephemeris and clock corrections for the corresponding GPS 
satellite when it is observed by at least four of the reference stations in 
25 accordance with: 

x k = (GvI G w )- 1 gJap" 

where: 

k represents the corresponding GPS satellite; 
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x* is a state vector containing the ephemeris and clock corrections for 
the corresponding GPS satellite; 

Ap k is a vector containing each synchronized pseudo-range residual for 
the corresponding GPS satellite, each synchronized pseudo-range residual in 
5 the vector Ap k corresponding to one of the reference stations that observes the 
corresponding GPS satellite; and 

Gyv is a matrix containing, for each reference station that observes the 
corresponding GPS satellite, a line of sight from it to the corresponding GPS 
satellite and an element that indicates that the corresponding synchronized 
10 pseudo-range residual in the vector Ap k should be adjusted by the clock 
correction for the corresponding GPS satellite. 

16. A GPS reference system as recited in claim 9 wherein: 

the network of reference stations comprises at least four reference 
15 stations; 

each reference station is further configured to compute a line of sight to 
each of the GPS satellites it observes; 

the master station is also in communication with the reference stations 
to receive the computed lines of sights; 
20 each ephemeris and clock correction estimator is configured to compute 

the ephemeris correction for the corresponding GPS satellite when it is 
observed by no more than three of the reference stations in accordance with: 

Ar* = Gvl (G w Gj) 1 A k 

and compute the ephemeris correction for the corresponding GPS satellite 
when it is observed by at least four of the reference stations in accordance 
25 with: 

Ar'MGwGwK 1 GwA k 

where: 

k represents the corresponding GPS satellite; 
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Ar* is a vector containing the ephemeris correction for the corresponding 
GPS satellite; 

A k is a vector containing, for each secondary reference station that 
observes the corresponding GPS satellite, the difference between the 
5 synchronized pseudo-range residual for the corresponding GPS satellite 
computed by it and the synchronized pseudo-range residual for the 
corresponding GPS satellite computed by the primary reference station; and 

Gyv is a matrix containing, for each secondary reference station that 
observes the corresponding GPS, the difference between a line of sight from 
10 it to the corresponding GPS satellite and a line sight from the primary reference 
station to the corresponding GPS satellite; and 

each ephemeris and clock correction estimator is further configured to 
compute the clock correction for the corresponding GPS satellite in response 
to the computed ephemeris correction for the corresponding GPS satellite. 

15 

17. A method of providing GPS satellite ephemeris and clock corrections to 
users utilizing a network of reference stations and a master station, the method 
comprising the steps of: 

with each reference station, receiving GPS signals from GPS satellites 
20 it observes and in response computing a pseudo-range residual for each 
observed GPS satellite; 

providing the pseudo-range residuals computed with each reference 
station to the master station; 

with the master station: 
25 in response to the provided pseudo-range residuals, computing 

clock differences between the reference stations and removing the clock 
differences from the provided pseudo-range residuals so as to synchronize 
them; 

for each GPS satellite observed by the reference stations, 
30 computing an ephemeris correction and a clock correction for the 
corresponding GPS satellite in response to the synchronized pseudo-range 
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residuals for the corresponding GPS satellite and independent of the 
synchronized pseudo-range residuals for the other GPS satellites observed by 
the reference stations; and 

transmitting the ephemeris and clock corrections for each GPS 
5 satellite observed by the reference stations. 

18. A method as recited in claim 17 wherein the ephemeris and clock 
correction computing step computes the ephemeris and clock corrections for 
the corresponding GPS satellite when the corresponding GPS satellite is 

10 observed by at least one of the reference stations. 

19. A method as recited in claim 18 wherein the network of reference 
stations comprises no more than three reference stations. 

15 20. A method as recited in claim 18 wherein the network of reference 
stations comprises at least four reference stations. 

21. A method as recited in claim 17 wherein: 

the network of reference stations comprises no more than three 
20 reference stations; 

for each GPS satellite observed by the reference stations, the ephemeris 
and clock correction computing step computes the ephemeris and clock 
corrections for the corresponding GPS satellite in accordance with: 

X k = Gj (G w Gjr 1 Ap k 

where: 

25 k represents the corresponding GPS satellite; 

x* is a state vector containing the ephemeris and clock corrections for 
the corresponding GPS satellite; 

Ap k is a vector containing each synchronized pseudo-range residual for 
the corresponding GPS satellite, each synchronized pseudo-range residual in 
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the vector Lp k corresponding to one of the reference stations that observes the 
corresponding GPS satellite; and 

G w is a matrix containing, for each reference station that observes the 
corresponding GPS satellite, a line of sight from it to the corresponding GPS 
5 satellite and an element that indicates that the corresponding synchronized 
pseudo-range residual in the vector Ap k should be adjusted by the clock 
correction for the corresponding GPS satellite. 

22. A method as recited in claim 17 wherein: 
10 the network of reference stations comprises no more than three 

reference stations and includes a primary reference station and no more than 

two secondary reference stations; 

for each GPS satellite observed by the reference stations, the ephemeris 

and clock correction computing step computes the ephemeris correction for the 
15 corresponding GPS satellite when it is observed by the primary reference 

station and at least one of the secondary reference stations in accordance with: 

Ar k =Gv5(G w Gw)- 1 A." 

where: 

k represents the corresponding GPS satellite; 

Ar* is a vector containing the ephemeris correction for the corresponding 
20 GPS satellite; 

A K is a vector containing, for each secondary reference station that 
observes the corresponding GPS satellite, the difference between the 
synchronized pseudo-range residual for the corresponding GPS satellite 
computed by it and the synchronized pseudo-range residual for the 
25 corresponding GPS satellite computed by the primary reference station; 

Gw is a matrix containing, for each secondary reference station that 
observes the corresponding GPS, the difference between a line of sight from 
it to the corresponding GPS satellite and a line sight from the primary reference 
station to the corresponding GPS satellite; and 



36 



WO 96/22546 



PCT/US96/00672 



for each GPS satellite observed by the reference stations, the ephemeris 
and clock correction computing step computes the clock correction for the 
corresponding GPS satellite in response to the computed ephemeris correction 
for the corresponding GPS satellite. 

23. A method as recited in claim 17 wherein: 

the network of reference stations comprises at least four reference 
stations; 

for each GPS satellite observed by the reference stations, the ephemeris 
10 and clock correction computing step computes the ephemeris and clock 
corrections for the corresponding GPS satellite when it is observed by no more 
than three of the reference stations in accordance with: 

x k = Gj (G w Gw)- 1 Ap k 

and computes the ephemeris and clock corrections for the corresponding GPS 
satellite when it is observed by at least four of the reference stations in 
15 accordance with: 

x'MGvjGwr'GjAp* 

where: 

k represents the corresponding GPS satellite; 

x k is a state vector containing the ephemeris and clock corrections for 
the corresponding GPS satellite; 
20 Ap K is a vector containing each synchronized pseudo-range residual for 

the corresponding GPS satellite, each synchronized pseudo-range residual in 
the vector Ap k corresponding to one of the reference stations that observes the 
corresponding GPS satellite; and 

Gyy is a matrix containing, for each reference station that observes the 
25 corresponding GPS satellite, a line of sight from it to the corresponding GPS 
satellite and an element that indicates that the corresponding synchronized 
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pseudo-range residual in the vector Ap k should be adjusted by the clock 
correction for the corresponding GPS satellite. 

24. A method as recited in claim 17 wherein: 
5 the network of reference stations comprises at least four reference 

stations including a primary reference station and at least three secondary 
reference stations; 

for each GPS satellite observed by the reference stations, the ephemeris 
and clock correction computing step computes the ephemeris correction for the 
10 corresponding GPS satellite when it is observed by no more than three of the 
reference stations in accordance with: 

Ar k = Gvl (G w Gyl)-' A* 

and computes the ephemeris correction for the corresponding GPS satellite 
when it is observed by at least four of the reference stations in accordance 
with: 

Ar' = (GvjGv,)- 1 Gyl A k 

15 where: 

k represents the corresponding GPS satellite; 
Ar* is a vector containing the ephemeris correction for the corresponding 
GPS satellite; 

A k is a vector containing, for each secondary reference station that 
20 observes the corresponding GPS satellite, the difference between the 
synchronized pseudo-range residual for the corresponding GPS satellite 
computed by it and the synchronized pseudo-range residual for the 
corresponding GPS satellite computed by the primary reference station; 

Gw is a matrix containing, for each secondary reference station that 
25 observes the corresponding GPS, the difference between a line of sight from 
it to the corresponding GPS satellite and a line sight from the primary reference 
station to the corresponding GPS satellite; and 
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for each GPS satellite observed by the reference stations, the ephemeris 
and clock correction computing step computes the clock correction for the 
corresponding GPS satellite in response to the computed ephemeris correction 
for the corresponding GPS satellite. 
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